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Abstract 

Hole-assisted optical fibers are microstructured fibers which guide light using a conventional high-index core while 
using a small number of air holes in the cladding to achieve desirable dispersive and/or polarization properties of the 
fiber. In this paper, we analyze the characteristics of the hole-assisted fibers using the multipole method of White et al. 
[Opt. Lett. 26 (2001) 1660], which is a rigorous boundary value method. Such properties as dispersion and birefringence 
are numerically investigated for fibers with different air hole patterns and sizes. Hole-assisted optical fibers exhibit low 
transmission loss and easily tailorable dispersion and birefringence, and therefore, would be valuable for dispersion 
management applications both in linear (WDM) and nonlinear (soliton) fiber telecommunications. © 2002 Published 
by Elsevier Science B.V. 



1. Introduction 

Air-silica microstructured optical fibers 
(ASMFs) (also termed photonic crystal fibers or 
holey optical fibers) have attracted a great deal of 
recent interest [1-9]. ASMFs consist of an ar- 
rangement of air holes running through the entire 
length of the fibers. Compared to conventional 
optical fibers, ASMFs exhibit properties such as a 
wide range of single-mode operation, tailorable 
dispersion, and easily controllable effective mode 
area. ASMFs are ready to lend themselves for 
applications such as dispersion compensation [4- 
6], soliton formation [7,8] and visible continuum 
generation [9]. 
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However, the high transmission loss of ASMFs 
is one of the problems that hinder many practical 
applications. The transmission loss in the 1.55 urn 
wavelength band has been typically as high as a 
few tens of dB/km. Although progresses in fab- 
rication of low-loss holey fibers with losses of a 
few dB/km have been recently reported [10,11], 
the losses are still significantly higher than those 
of the conventional fibers (typically 0.2-0.5 dB/ 
km). In recent efforts to reduce the transmission 
loss of ASMFs, Hasegawa et al. [12-14] proposed 
and fabricated "hole-assisted lightguide fibers" 
with losses comparable to conventional fibers. 
The low-loss fiber has a high index core, a low 
index cladding, and a few air holes surrounding 
the core as shown in Fig. 1. It can be treated as a 
conventional fiber with a few air holes intro- 
duced into the cladding, therefore, a variation 
lying between the conventional fiber and the ho- 
ley fiber. 
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Fig. 1. Schematic of the air-hole-assisted optical fiber. The 
high-index core of radius R 0 is surrounded by N(N = 6 shown) 
circular air holes. 



The optical properties of ASMFs have been 
modeled by a number of methods including both 
approximate methods such as the effective index 
method [1,15] and rigorous methods such as 
plane-wave-method [16,17], localized-function 
method [18,19], finite element method [20]. All 
these rigorous approaches start from the trans- 
verse field component and solve for propagation 
constant /? and associated modal fields. For ex- 
ample, the localized-function method expands the 
transverse electric or magnetic field in terms of a 
Hermite-Gauss basis and solves the resultant al- 
gebraic eigenvalue equation for /?. Recently, 
multipole method has been used in modeling the 
confinement loss of ASMFs ([22] and references 
therein), in which the EM properties of ASMFs 
are obtained by starting from the longitudinal 
field components and solving a boundary value 
problem. Hasegawa et al. [13] have used the full- 
vector finite element method to model the air- 
assisted fiber. In this paper, we use the multipole 
method to investigate the optical properties of the 
air-hole-assisted fiber. 

In Section 2, the formation of the multi- 
pole method is reviewed. In Section 3, the ef- 
fective index, group dispersion and birefringence 
of the fundamental mode are presented for fibers 



with different air hole number and size. Com- 
parisons are also made between the air-assisted 
fiber and the conventional fiber and between the 
multipole method and the localized-function 
method. 



2. Formulation 

The multipole method starts from the longitu- 
dinal field components E z and H 2 instead of the 
transverse field components, just like in the case of 
conventional optical fibers. The longitudinal field 
satisfies the Helmholtz equation 

v;Q) +( w-rt(J)-o, (.) 

where ft is the propagation constant in z direction, 
and n = n{x,y) is the transverse refractive index 
profile of the fiber. The longitudinal field compo- 
nents determine the transverse components from 
Maxwell's equations: 

E ' = F^?( v ' E '- z "? ixv,H ')' p) 
h '=A( v ' h ' + ^ 2,<v ' e -)' <3) 

where Zo = is the impedance of free space. 

From the boundary conditions at all the inter- 
faces (i.e., the continuity of tangential components 
of electric and magnetic fields), we can solve for 
the propagation constant jS and corresponding 
fields. 

We consider the fiber with a high-index guiding 
core and N circular air holes as illustrated in Fig. 
1. The high-index core has refractive index n 0 and 
radius Ro, The y'th air hole has index of unity and 
radius Rj. The cladding has refractive index 
n c (n c < n Q ) and it is assumed to be infinite.Within 
each circular region, the E Zi H z are expressed in the 
form 

E*J = £ exp(im0,), (4a) 

m~-oo 

H 2J = Z 0 -' £ Wipjrj) exp(im0,), (4b) 

m=-oo 
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where 



A.(o>/?-^) 



otherwise, 



(5) 



and (r,, 0 7 ) is the local polar coordinate with the 
yth circular region. 

In the cladding, the longitudinal field compo- 
nents have the form 



/=0 m=-oo 

//« = Z 0 "' £ £ DgKmhcTi) exp(im0,), 



(6a) 
(6b) 



/-0 m— -oo 



where <? c = y^/? 2 - k 2 n\. 

In the above expressions, the multiple coordi- 
nate systems associated with each circular region 
are used. Using the addition theorem for Bessel 
function K m [22], we can rewrite the E x and H ZQ in 
the yth coordinate system by 



m=-oo 

+ C^ m ( 9c o)]exp(if«0;) ) 

oo 
m=-oo 

+ exp(im0,), 

with 

oo 

¥7 9=-<» 



(7a) 



(7b) 



exp 



9 =-oo 



exp 



(8a) 



(8b) 



where (rj,0j) is the polar coordinate of the /th 
origin in the y'th polar coordinate system. 



From the longitudinal field components, we can 
get the transverse tangential components using 
Eqs. (2) and (3) 



E Bj = 



E 



-kPjB^ipjrj) 



exp(im9j), 



(9a) 



iz 0 -' 



— — ■ — «• Y kn)pjA^F^(p jrj ) 
k 2 n] - p 2 [ jFj m m W Jl 

+ ^(im)B^F^{p J rj)\ e xp(ime J ), 



+ C^ M ( 9c i V )]-A 9c [KW/:( gc0 ) 
+ £>W^(?cr / )]}exp(im^) ) 

• 7 -l OO ( 

"* = FTT^ £ \kn 2 MU^I' m (q c rj) 

K U c P m=-oo ^ 

+ C^K(g c rj)]+^{im)[^I m (q e rj) 
+ DMK m (q c rj)]}exp(imdj). 



(9b) 



(9c) 



(9d) 



Applying Eqs. (4a)-(9d) to the boundary con- 
ditions at all the N + 1 circular interfaces, we ob- 
tain an infinite system of linear equation in terms 
of the unknown coefficients C$ and D$(j = 
0, 1, . . . ,JV). I* can be written in a matrix form as 
MC = 0, (10) 

where C is the vector composing of unknown 
constants Cjp and D$, and M is the coefficient 
matrix whose determinant vanishes for a nontriv- 
ial solution of C. In order to numerically calculate 
the propagation constant P of the guided modes, 
we truncate M into a finite square matrix. If we 
limit m in the range from —L to L, then the di- 
mension of the resultant matrix is 2(N+1) 
(21 + 1) x 2(Af + 1)(2£ + 1). The propagation 
constants of guided modes are then obtained by 
iteratively scanning for dips of the determinant of 
the truncated matrix M while fi/k is varied from to 
n c to no. 
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Fig. 2. Calculated (a) effective index and (b) group velocity dispersion (GVD) for air-assisted fibers with different number of identical 
air holes symmetrically surrounding the core. The core is doped with 10% mol Ge0 2 and has a radius R 0 of 2 urn. The core-to-hole 
spacing A is 5 um. Air holes have a radius R a of 2 urn. 



3. Results and discussion 

We first consider the microstructured optical 
fiber as illustrated in Fig. 1. The central core is 
Ge0 2 doped silica and is symmetrically sur- 
rounded by identical air holes in the pure silica 
cladding. It has fewer air holes than other ASMFs 
since the higher index core can efficiently confine 
light and thus enables low-loss propagation. When 
the air holes are small or far from the core, the 
fiber behaves like a conventional fiber. However, 
the air holes will significantly alter the dispersion 
of the fiber, especially when they are large or close 
to the core. 

In Fig. 2, we show the calculated propagation 
constants and group velocity dispersions (GVD) of 
the fundamental guided modes in the fibers with 3, 

4, 6, and 8 identical air holes. The corresponding 
conventional fiber (N = 0) is also shown for 
comparison purpose. It can be seen that the pres- 
ence of air holes shifts the zero dispersion wave- 
length A D to shorter wavelength and that more air 
holes can increase anomalous dispersion in the 
long wavelength region. Alternatively, large air 
hole or being closer to the core can also increase 
anomalous dispersion as shown in Fig. 3. 

The modal fields, Poynting vector component 
of the fundamental mode (x-polarized) are shown 



in Fig. 4 for a 6-air-hole fiber at wavelength of 
1 .5 n.m. As can be seen, the air holes do not change 
the energy distribution very much; average energy 
flow density (S z ) is essentially the same as a con- 
ventional fiber. This feature should make the fiber 
easy to splice with conventional single mode fiber 
with little loss. However, both the propagation 
constant and GVD are changed significantly by 
the air holes. 




Wavelength (um) 



Fig. 3. Calculated GVD for 6-air-hole fibers with different air 
hole sizes. The core is doped with 10% mol Ge0 2 and has a 
radius R Q of 2 um. The core-to-hole spacing A is 5 um. 
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Fig. 4. (a) \E z \ t (b) \H t \, (c) E„ and (d) (S z ) for the fundamental mode (^-polarized) of a 6-air-hole fiber at wavelength of 1.5 urn (core 
radius R 0 = 2 urn; 10% mol Ge0 2 doping in the core; air hole radius R a = 2 urn; hole-hole spacing A = 5 urn). 



The fiber can be made to be polarization 
maintaining by inducing high birefringence 
through changing either the position or size of 
some of the air holes. As an example, here we 
consider changing the size of some air holes as il- 
lustrated in Fig. 5(a). In Fig. 5(b), we show the 
calculated birefringence B when the radius R b of 
the two opposite air holes is varied. 

We notice that, for long wavelengths, the ef- 
fective index of the fundamental mode cannot be 
found via the searching technique mentioned in 
Section 2, because the modal fields extend more 
significantly into the air hole the effective index of 



the mode f$/k is smaller than n c . In this limit, /? is 
complex and the guided mode becomes lossy. 
However, the multipole method is still capable of 
solving the problem by using more general for- 
mulations and proper boundary conditions as in 
[21]. As an alternative, we use the localized-func- 
tion method [18,19] for the long wavelengths. We 
show the results from both methods in Fig. 6. They 
agree very well, especially in the short wavelength 
region. 

The multipole method verifies the degeneracy of 
the fundamental mode of the air-hole-assisted fiber 
with Af-fold (N > 3) symmetry (More general 
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Fig. 5. Birefringence of a 6-air-hole fiber when the size of the 
two opposite air holes varies, (a) Schematic of the fiber con- 
sidered in the calculation. R a is fixed, and R b varies, (b) Cal- 
culated birefringence for fibers shown in (a) with different R b . 
The core is doped with 10% mol Ge0 2 and has a radius R 0 of 
2 jim. R a - 2 um. The core-to-hole spacing A is 5 um. 



discussions of the symmetry and degeneracy of 
microstructure optical fibers have been made in 
[23,24]). For a 3-fold symmetric hole-assisted fiber, 
in Fig. 7 we show the determinant of the matrix M 
when the effective index fi/k is scanned from n c to 
n 0 . The presence of only one sharp dip indicates no 
splitting of the fundamental mode and verifies that 
this method introduces no numerical breaking of 
the degeneracy. 
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Fig. 6. Comparison between the multipole method and the 
localized-function method for a fiber with six symmetrically 
arranged identical air holes. The core is doped with 5% mol 
Ge0 2 and has a radius R 0 of 3 um. The core-to-hole spacing A 
is 6 um. Air holes have a radium R a of 2 um. 
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Fig. 7. The determinant of M vs. effective index for a fiber 
with three symmetrically arranged identical air holes. The inset 
shows details around the dip. The core is doped with 5% mol 
Ge0 2 and has a radius R 0 of 3 um. The core-to-hole spacing A 
is 6 um. Air holes have a radius R a of 2 um. 



4. Conclusion 

We have analyzed the guiding properties of a 
special type air-silica microstructured fiber that 
has a high-index guiding core and a single ring of 
air holes in the cladding, by solving a rigorous 
boundary value problem. The group velocity 
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dispersions of the fiber with different number of air 
holes and different air hole sizes are numerically 
obtained. The birefringence of the fiber by 
changing the air hole size is also investigated. With 
important features such as low transmission loss, 
tailorable dispersion and polarization, the air- 
hole-assisted optical fiber may be of great promise 
in lightwave applications. 
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